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ABSTRACT 

In a recent preprint, iHearin et al.l ( 120121 HI 2) suggest that the halo mass-richness calibration of clusters 
can be improved by using the difference in the magnitude of the brightest and the second brightest galaxy 
(magnitude gap) as an ad ditional observable. They claim that their results are at odds with the results from 
Paraniape & Sheth (120121 PS 12) who show that the magnitude distribution of the brightest and second brightest 
galaxies can be explained based on order statistics of luminosities randomly sampled from the total galaxy 
luminosity function. We find that a conditional luminosity function (CLF) for galaxies which varies with halo 
mass, in a manner which is consistent with existing observations, naturally leads to a magnitude gap distribution 
which changes as a function of halo mass at fixed richness, in qualitative agreement with H12. We show that, 
in general, the luminosity distribution of the brightest and the second brightest galaxy depends upon whether 
the luminosities of galaxies are drawn from the CLF or the global luminosity function. However, we also show 
that the difference between the two cases is small enough to evade detection in the small sample investigated 
by PS 12. This shows that the luminosity distribution is not the appropriate statistic to distinguish between 
the two cases, given the small sample size. We argue in favor of the CLF (and therefore HI 2) based upon its 
consistency with other independent observations, such as the kinematics of satellite galaxies, the abundance 
and clustering of galaxies, and the galaxy-galaxy lensing signal from the Sloan Digital Sky Survey. 

Subject headings: cosmology: theory ■ 
halos 



methods: numerical - galaxiesxlusters - galaxies:evolution - galaxies: 



1. INTRODUCTION 

The abundance of halos at the massive end is sensitive to 
variations in cosmological parameters such as the matter den- 
sity parameter (£2 m ) and the amplitude of the power spectrum 
of density fluctuations in the Universe (characterized by erg). 
Therefore observations of the abundance of galaxy clusters 
which reside in such massive halos ca n be used to constrain 
these cosmological parameters (see e.g lVikhlinin et aDl2009t 
iMantz et all 120 lot iRozo et al.lfeOlOl) . Additionally, such ob- 
servations can also be used to constrain the phenomenological 
behavior of dark energy and test the nature of gravity through 
measurements of the growth of structure, questions which are 
of fundamental importan ce for cosmologists today (see e.g., 
iRapetti et al.1 120 lOT 1201 21). Photometric surveys such as the 
Dark Energy Survev Ubark energy survey collaborationl2005 , 
DES), the Hyper-Suprime Cam survey dMivazaki et al.ll2006 , 
HSC) in the immediate future and the Large Scale Synoptic 
Telescope Survey (LSST) in the near future, will result in a 
large catalog of optically-selected galaxy clusters, which can 
be used to answer these important questions. 

Identifying the galaxy cluster observables in optical surveys 
which tightly correlate with halo mass, establishing the scal- 
ing relations between these observables and halo mass and 
the scatter in these relations are all crucial steps in order to 
achieve the scientific goals. It is well known that the num- 
ber of cluster mem bers (also called richness) correlates with 
halo mass (see e.g, iBecker et alj|2007l Uohnston et al.ll2007t 
iSheldon et al.l 120091) . and so does the lu minosity (or stellar 
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tent in the group (see e.g.. lYang et a l. 2007, for results based 
on abundance matching). However, these scaling relations 
have considerable scatter, and therefore combining multiple 
observables, especially those which come without additional 
observation al costs, is an impor tant task. 

Recently. IHearin et alj ([20 12) suggested that at fixed rich- 
ness, the magnitude difference (also called magnitude gap 
or equivalently the luminosity ratio) between the brightest 
and the second brightest galaxy, contains information about 
halo mass. By using a simple subhalo abundance match- 
ing prescription they showed that at fixed richness, the small 
(large) magnitude gap systems are expected to preferentially 
reside in less (more) massive halos. To test the proposition 
with real data, the a uthors used the galaxy group catalog of 
iBerlind et al.l (|2006) and showed that at fixed velocity dis- 
persion (proxy for halo mass), the average richness of small 
magnitude gap systems is significantly larger than the average 
richness of the large magnitude gap systems, and the differ- 
ence is larger than that expected if the luminosities of galaxies 
in every group were drawn randomly from the galaxy lumi- 
nosity function. 

Theoretically, it is expected that dynamical friction causes 
brighter satellites in massive halos to merge with the central 
galaxy, increasing the magnitude gap between the brightest 
satellite and the central galaxy. This suggests that central 
galaxies are expected to be special, they occupy the deepest 
portion of the potential well, where they grow by feeding on 
the satellites that are dragged to the center of halos by dynam- 
ical friction. Whether this is indeed the case, or whether the 
luminosity of the central galaxies is just a matter of chance, is 
a ma tter which can be se tt led wi th observations. 

In iParaniape & She till J2012I) . the authors examine this 
question by looking at the luminosity distribution of the 
brightest and the second brightest galaxy and t he magnitude 
gap between the two using the group catalog of IBerlind et alj 
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( 120061) . They find that the luminosity distributions are con- 
sistent with the distribution of the brightest and the second 
brightest of N random draws from the galaxy luminosity func- 
tion, where N is the richness of a given group. On face value, 
this would imply that there is nothing special about the bright- 
est galaxy in a given group, and that it is just a matter of 
chance that any galaxy becomes the brightest in a given group. 
These results imply that the magnitude gap should not con- 
tain any more information about the halo mass, than that con- 
taine d in the richness, i n apparent contrast with the results 
from lHearin et al.1 (120121) . which are based on the same group 
catalog. 

In this paper, we attempt to clarify this issue, by predict- 
ing the magnitude gap based upon the conditional luminosity 
function (CLF), which describes the halo occupation distribu- 
tion of galaxies in a halo of given mass. The CLF and its vari- 
ation with halo mass has been calibrated using a wide variety 
of observations such as the abundance of galaxies, their clus- 
tering and the galaxy-gala xy lensing signal measured from the 
Sloan digital sky survey (I York et al.ll2000l SDSS). We show 
that if galaxies occupy halos according to the CLF, it is nat- 
ural to expect that the magnitude gap depends upon the halo 
mass at fixed richness. We also show that the luminosity dis- 
tributions of the brightest and the second brightest galaxies 
are predicted to be different from those obtained by random 
draws from the luminosity function. However, detecting this 
small difference just using the luminosity distribution will 
require sample sizes whi ch are larger than the one used by 
IParanjape & Shethl (120121) . 

We also note that in their paper, IParanjape & Shethl (120121) 
investigate the luminosity-weighted marked correlation func- 
tion and show that its radial dependence implies that the lu- 
minosities of the brightest galaxies are not a matter a chance. 
They conclude that their results falsify the hypothesis that the 
luminosities of the brightest and the second brightest galaxy 
are drawn from the global luminosity function (i.e., without 
any dependence on halo mass or environment) and that the 
luminosity distribution alone is not an appropriate discrimi- 
nant to investigate this issue. In this paper, our results based 
on the CLF will strengthen this argument. 

This paper is organized as follows. In Section [2] we de- 
scribe the CLF framework and give analytical expressions for 
the magnitude gap distribution based upon the CLF. In Sec- 
tion [3j we show the magnitude gap distribution from Monte 
Carlo simulations based on the CLF and compare the re- 
sults to the analytical expression presented in Section [2] We 
also investigate the dependence of the magnitude gap upon 
the richness in a group and the assumed CLF parameterisa- 
tion. In Section |U we construct mock galaxy catalogs based 
upon galaxy luminosities sampled from (a) the CLF and (b) 
the overall galaxy luminosity function and compare the lumi- 
nosity distributions of the brightest and the second brightest 
galaxy in these two catalogs. Finally, we summarize our re- 
sults in Section|5] 

For the purposes of this paper, we adopt the following con- 
vention. We refer to galaxies as centrals (satellites), if they 
are drawn from the CLF which is specific to the central (satel- 
lite) galaxies (see Eqs.|2]and|4]i. As our fiducial model, we 
assume that central galaxies are also the brightest galaxies in 
the halo. Therefore, in the fiducial case, the magnitude gap is 
the difference in magnitudes between the central galaxy and 
the brightest satellite. However, we will also investigate cases, 
when the sa tellites are allowed to be brighter than the central 
galaxy (see ISkibba et al.ll2.Ql ll for observational evidence of 



such a possibility). 

2. CONDITIONAL LUMINOSITY FUNCTION 

The conditional luminosity function, denoted by <1>(L|M), is 
defined to be the average number of galax ies of luminositie s 
L ± dL/2 that reside in a halo of mass M (lYang et al.ll2003l) . 
The average number of galaxies in a given halo of mass M can 
be found by simply integrating the CLF over the luminosities 
of interest, e.g., the average number of galaxies with lumi- 
nosities between L m j n and L max that reside in a halo of mass 
M is given by 



(N)m(Lowi> ^max) 



r 



<D(L|M)dL . 



(1) 



For convenience, the CLF is divided in to a central galaxy 
component (<i> c [L|M]) and a satellite galaxy component 
(0> s [L|M]). 

We assume that the distribution O c (L|M) is described by 
a lognormal distribution with a scatter, <x c , that is indepen- 
dent of halo mass, co nsistent with the fi nd ings f rom studies 
of satellite kine matics (iMore et al Jl2009alfbL f20 1 ll) and galaxy 
group catalogs (lYang et alJ 120091) . 
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The dependence of the logarithmic mean luminosity, logL c , 
on halo mass is given by 



log L C (M) = log 



Lo 



(M/MiY 



[1 + (M/M0] 



7l -72 



(3) 



Four parameters are required to describe this dependence; two 
normalization parameters, Lq and M\ and two parameters j\ 
and 72 that describe the slope of the L C (M) relation at the low 
mass end and the high mass end, respectively. 

The satellite CLF, O s (L|M) is assumed to be a Schechter- 
like function, 



<D s (L|M)dL = O* 



/. 



exp 



Lj 



AL 



(4) 



Here L t (M) determines the knee of the satellite CLF and is 
assumed to be a factor / s times fainter tha n L C (M). Motivate d 
by results from the SDSS gr oup catalog of 
we set f s = 0.562 (see also iReddick et al] 



Yang et alJ 
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!2008) 
2, and 

assume that the faint-end slope of the satellite CLF is indepen- 
dent of halo mass. The logarithm of the normalization, O* is 
assumed to have a quadratic dependence on log M described 
by three free parameters, bo, b\ and bi\ 



log<D* = bo + bi (logM- 12) + MlogM- 12) 2 



(5) 



Note that this functional form does not have a physical moti- 
vation; it merely provides an adequ ate description of the re- 
sults obtained by lYang et"ail (120081) from the SDSS galaxy 
group catalog. The parameters of the conditional luminos- 
ity function and their variation with halo mass can be con- 
strained by using observations of the abundance, the cluster- 
ing and the galaxy-galax y lensing signal me a sured from the 
Sloan Digital Sky Survey ( IMore et al.l2012allHlc"acciato et ail 
120121) . In what follows,' we will use the following values for 
the CLF parameters: L = 10 995 /T 2 L , M\ = lO 11 ' 27 /^ 1 M , 
cr c = 0.156, 7i = 2.94, 72 = 0.244, a s = -1.17, b = -1.42, 
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Fig. 1 . — The distribution of the difference in magnitudes between the brightest and the second brightest galaxy predicted by simulations in which we populate 
galaxies in halos of different mass according to the conditional luminosity function. The solid histograms show the distribution of magnitude gaps in halos of 
different mass (shown using different colors) for our fiducial model in which we assume that the central galaxy (defined to be drawn from the central CLF) is the 
brightest in the halo. The solid curves show the analytical prediction based on Eq.|6] The dashed histograms show the corresponding result for the case when we 
allow satellites to be brighter than the central galaxy. 

b\ = 1.8 2, and b 2 - -0.30, c onsistent with the results pre- 
sented in lCacciato et al.l ( 120121) . 

If the luminosities of galaxies in a halo are drawn in an 
uncorrelated fashion, the probability that a halo of mass M 
and richness N has a magnitude gap, Am, or equivalently the 
luminosity ratio, fi, between the brightest satellite galaxy and 
the central galaxy in a halo of mass M is then given b>0 

dL' P s (L'\M)P c (L'/f L \M) 



function, respectively, such that 
(N c ) M (LuL 2 ) = l - erfc 



-erfc 



(N s ) M (LuL2) = — II 
P 



log Li - logL c 

V2<T C 

logL 2 -logL, 

V2cr c 
I 



a s + l (L-V 



x[P s (<L'\M)T 



(6) 



Here, the probabilities P X (L'\M) and P x (< L'\M) are defined 
such that 



P \L 
a s + l (L 2 \ p 
L, 



(9) 



(10) 



P,{L'\M) = 
PA< L'\M) = 



O x (L'|M) 



(^Vx)m (^min; ^max) 



(7) 
(8) 



The probability of a halo to have a certain mass, given its 
richness and the magnitude gap can be obtained from Eq. [6] 
and the Bayes' theorem, 



P(M\f L ,N) 



P(f L \M,N)P(M\N) 
P(/l\N) 



(11) 



where the symbol x can either stand for central (c) or satellite 
(s). The quantities (N x )m in the relevant luminosity intervals 
can be obtained by replacing 0(L|M) by <f> x (L\M) inside the 
integral in Eq. ([TJ- For central galaxies we choose L max = oo. 
In our model, we assume that the central galaxies are always 
the brightest in the halo. Therefore, in the case of satellites, 
we use the luminosity of the central under consideration as 
the upper limit, i.e., L max = 

The integrals for (N c )m and (N s )m can be written in terms of 
the complementary error function and the incomplete gamma 



and as expected it depends upon the mass-richness relation 
via the probability distribution P(M\N). The probability dis- 
tribution within a given bin of richness [Ni , N 2 ] is given by 



P{M\f L ,Nx <N<N 2 )= J) 



P(f L \M, N)P(N\M)P(M) 



N=N, 



P(/l\N) 



4 Note that our expression differs from Paranjape & Sheth 12012) because 
in our case the central galaxy luminosity is assumed to be sampled from a 
probability distribution which differs from the distribution from which the 
satellites are sampled from. 



(12) 

Finally, the distribution of the magnitude gap at fixed halo 
mass (without regard to the richness) is given by 

DO 

P(fi\M) = ^ P(f L \N, M)P(N\M) . (13) 

In what follows, we will also investigate the effect of allow- 
ing satellite galaxies to be brighter than the central galaxies 
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in their halo. The analytical expressions for predicting the 
magnitude gap distribution in this case are presented in the 
appendix. 

3. RESULTS FROM SIMULATED SAMPLE 

We now demonstrate explicitly that for fixed richness, the 
CLF predicts that the magnitude gap in a given group of 
galaxies depends upon the mass of the halo in which these 
galaxies reside. The CLF varies with halo mass and there- 
fore it is not surprising that this indeed is the case. For this 
purpose, we generate Monte-Carlo samples of galaxies that 
populate halos according to the CLF in the following manner. 

For a halo of given mass, we first draw the luminosity of 
its central galaxy from (S> cen (L\M), given by Eq. (|2). In order 
to avoid the existing correlation between halo mass and rich- 
ness affecting our conclusions, we fix the number of satellites 
A^ S at = 20. For each of the iV sat satellites, we then draw a lu- 
minosity from the satellite CLF Q sat (L\M), given by Eq. ©. 
While drawing the satellite luminosities, we adopt a lumi- 
nosity threshold, L m ; n , corresponding to 01 M r — 51og/2 = 
-19 (here 0A M r indicates the SDSS r-band magnitude, K- 
corrected to z = 0.1; see Blanton et al. 2003). As mentioned 
before, we also assume that the satellites are always fainter 
than the central galaxy drawn for a given halo. 

The resultant distribution of the magnitude gaps is shown 
in Fig. [T] using a solid histogram for a wide range of halo 
masses. It can be clearly seen that the distribution of the mag- 
nitude gaps depends upon halo mass. We use Eq.[6]to predict 
this distribution analytically and compare it to the results from 
our simulations. The result of this analytical calculation are 
shown as solid curves in Fig. Q] which agrees well with the 
magnitude gap distribution from our simulations. 

For low mass halos, the distribution of magnitude gaps is 
peaked at zero. However, this peak shifts away from zero as 
we move to larger halo masses. This figure establishes that if 
galaxies populate halos according to the conditional luminos- 
ity function (which is supported by several observations such 
as galaxy group catalog, and the observations of abundance 
clustering and galaxy galaxy lensing from SDSS), the magni- 
tude gap should have more information about the halo mass, 
in addition to that conveyed by richness alone. At fixed rich- 
ness, higher mass halos tend to have la rger magnitude gaps, 
in agreement with iHearin et al.l (120121) . Our result that the 
magnitude gap distribution for low mass halos is peaked at 
zero, and shifts to larger magnitude gaps for larger mass ha- 
lo s, may appear to be exactly opposite of the result presented 
in Ivan den Bosch et al.l (120071) . However, note that the mag- 
nitude gap distributions we present are at fixed richness and 
halo mass (P[fi\N, M]), while the magnit ude gap distributions 
shown in the different panels in fig. 5 of Ivan den Bosch et al.l 
( 120071) correspond to groups with varying richness (thus cor- 
responding to P{fi\M), see Eq. [13) , due to the underlying 
mass-richness relation. We will shortly consider the effect of 
changing richness on the magnitude gap distribution. 

First, we consider the effect of relaxing the assumption that 
the centrals are the brightest in their halos. Note that in this 
case, the magnitude gap could be either between the cen- 
tral and the brightest satellite, or between the two brightest 
satellites, in case the halo has two or more satellites brighter 
than the putative central galaxy (see Appendix). The resul- 
tant magnitude gap is shown with a dashed histogram. For 
low mass halos it can be hardly distinguished from the case 
when we demand the central to be the brightest. It can be 
also seen that for all halo masses the distribution of magni- 



tude gaps for Ami2 > 0.5 is consistent with the case when 
the central galaxy is assumed to be the brightest. This is also 
expected since the satellite conditional luminosity function in 
our model dies exponentially at the bright end. Therefore if 
there is a satellite galaxy brighter than the central galaxy, the 
magnitude gap is not expected to be extremely large. There- 
fore, the few cases when the satellite galaxy is brighter, cause 
a small but noticeable increase in the probability distribution 
at the small magnitude gap end. 

We show the results of varying the number of satellites in 
Fig. [2] The left hand panel shows the magnitude gap distri- 
bution in halos of different mass, when the number of satel- 
lites equals 5, while the right hand panel shows the same for 
number of satellites equal to 64. As the number of satellites 
increases (decreases) the magnitude gap tends to be smaller 
(larger), as expected (and q ualitatively consistent with the re- 
sults presented in fig. 5 of Ivan den Bosch et al.l 12007). The 
analytical expectation (from Eq. [6) is shown as a solid curve; 
it describes the simulation results accurately, and is shown as 
a sanity check. 

The parameter p governs the exponential cut-off at the 
bright end of the satellite conditional luminosity function (see 
Eq.@). Based upon the analysis of offsets of the line-of-sight 
velocities and projected position of the brighte st galaxy rel- 
ative t o the mean of the other group members, Ski bba et al.l 
d201 ll) concluded that the value of p ought to be closer to 
unit y instead of the fidu cial value of 2 that we assume (see 
also Reddick et al. 2012). Therefore, we also show the effect 
of varying the parameter p on the magnitude gap distribution 
in Fig. [3] We have verified that the predictions based upon 
Eq. [6|that we show in the figure also agree with detailed sim- 
ulations. As expected, decreasing the value of p causes the 
satellite conditional luminosity function to fall less rapidly at 
the bright end which results in smaller magnitude gaps. 

Regardless of these details, it is clear that the results from 
this section establish that if galaxies populate halos according 
to the CLF, then at fixed richness the magnitude gap distribu- 
tion should depend upon the halo mass, in a ma nner which is 
qualitatively consistent with lHearin et al.l ( 120121) . 

4. LUMINOSITY DISTRIBUTION OF THE BRIGHTEST AND 
SECOND BRIGHTEST GALAXY 

We would like to now investigate the result presented in 
iParanjape & Shethl (120121) . They demonstrate that the lumi- 
nosity distribution of the brig htest and the second brightest 
galaxies in the group catalog o flBerlind et al.l (120061) is consis- 
tent with their expected distribution if the luminosity of galax- 
ies in each of the groups were randomly sampled from the 
global luminosity function of galaxies. To verify their result, 
we construct Monte-Carlo galaxy catalogs in which galaxy 
luminosities are drawn either from the conditional luminosity 
function or the overall luminosity function. 

We assume a standard flat ACDM cosmological model with 
matter density Q. m = 0.27, baryon density = 0.0469, 
Hubble parameter h = 0.7, spectral index n s = 0.95, and 
a matter power spectrum normalization of erg = 0.82. We 
sample a large number of haloes with masses M > 2.7 x 
10 13 r'M from the halo mass function expected for such 
co smology using the h alo mass function calibration presented 
by iTinker et "all (120081) . 

We construct a mock galaxy catalog (Catalog A) by pop- 
ulating the dark matter halos with model galaxies using the 
CLF with parameters described in $2] For each halo, we 
first draw the luminosity of its central galaxy from <& cen (L\M), 
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Fig. 2. — Distribution of magnitude gaps (similar to Fig.[T} for the case when the number of satellites is equal to 5 and 64 is shown in the left and right hand 
panel, respectively, and assuming that the central galaxy is the brightest in the halo. 



E 
< 



2 - 



1 ^ 



i — ■ 1 ■ 1 1 1 1 1 1 1 1 1 


i 1 i i i i 


\ N s =20 logM(h" 


'M )=13 - 


~\ logM(h" 


'M )= 15.25 


P= 


= 1 


P= 


=2,Fiducial - 


-. p = 


=3 






— ? V C~*<>. 

i i i i 1 i i i 









0.5 



1.5 



Am ia m sec brighl ni bright 

Fig. 3. — The dependence of the distribution of magnitude gaps on the 
parameter p which governs the exponential cutoff at the bright end in the 
conditional luminosity function of satellite galaxies, based on Eq. [6] The 
solid line corresponds to the fiducial case p = 2 while the dotted and the dot- 
dashed lines correspond to p = 1 and p = 3, respectively. The results for the 
two different halo masses are shown by different colors. 



given by Eq. (f2]). Next, we draw the number of satellite galax- 
ies, under the assumption that P(N sal \M) follows a Poisson 
distribution with mean given by Eq. (Q~|) with O replaced by 
<5 S , and we adopt a luminosity threshold, L m i n , corresponding 
to 01 M r -5 log/; = -20, similar to the threshold adopted by 
iParaniape & Shethl d2012l) . For each of the N sdt satellites in the 
halo of question, we then draw a luminosity from the satellite 
CLF <$) sat (L\M), given by Eq. (HJi and maintain the fiducial as- 
sumption that all satellites are fainter than the central galaxy. 
We restrict ourselves to halos with richness N > 12, which 
gives us a sample of 319482 halos. 

We construct an alternate catalog of galaxies (Catalog B) 
where the luminosities of member galaxies in each halo are 
drawn from the global luminosity function, <£>(L), 



<J>(L) 



-/ 



<&(L\M)n(M)dM , 



(14) 



where n(M) is the halo mass function. In practice, we ran- 
domly sample (with replacement) from the luminosities of 
galaxies in the entire previous catalog, while maintaining the 
richness of the group they belong to, thus effectively sampling 
the galaxy luminosities in every group from the global lumi- 
nosity function of galaxies. 

The luminosity distribution of the brightest and the second 
brightest galaxies in each halo for both the catalogs are shown 
in the upper left and right hand panels of Fig. [4] respectively. 
The peak of the magnitude distribution of the brightest galax- 
ies in Catalog A have a distribution which peaks at a slightly 
higher value of luminosity compared to Catalog B. On the 
other hand, the magnitude distribution of the second brightest 
galaxies shows a tail towards larger luminosities in Catalog B 
compared to that in Catalog A. However, the plot also shows 
that the differences are not that huge, and detecting such dif- 
ferences in the magnitude distributions will require a large 
sample of groups. 

From our large sample of Monte-Carlo groups, we now re- 
strict ourselves to selec ting sample sizes (~ 350) which are 
similar to those used bv lParanjape & Shethl (120121) . We show 
the results of one of the random realizations in the bottom 
panel of Fig. |4] We also obtain the corresponding cumulative 
distributions and use the Kolmogorov-Smirnov (KS) statis- 
tic to compare the distributions from the two catalogs. The 
p-values from the KS-test are indicated in the corresponding 
panels and these values imply that the luminosity distributions 
from the tw o catalogs, when downsam pled to the size of the 
catalog that IParaniape & Sheth| d2012l) use are consistent with 
each other. To show that this particular random realization 
is not a statistical fluke, in Fig. [5] we show the distribution 
of p-values from KS-tests carried out on 1000 random sam- 
ples si milar in size to the catalog used bv lParanjape & Shethl 
(120121) . The distribution of p-values from the KS-test peak 
at values larger than 0.1, which highlights the difficulty in 
distinguishing between the magnitude distributions from the 
two catalogs with a sma ll sample size. Th i s sugg ests that 
the group catalog used by lParanjape & Shethl (120 121) does not 
have enough number statistics, to detect the difference be- 
tween the luminosity distributions of the brightest (or the sec- 
ond brightest) galaxies in the cases corresponding to the two 
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Fig. 4. — Comparison between the luminosity distribution of the brightest and the second brightest galaxy in the halos present in the two mock galaxy catalogs 
are shown in the left and right hand side panels, respectively. Upper panels: The solid line shows the luminosity distribution when galaxies are populated in 
halos according to random draws from the CLF (Catalog A), while the dashed histogram shows the distribution when galaxies are populated according to random 
draws from the global luminosity fu nction (Catalog B), mainta ining the richness of halos. Bottom panels: Same as the upper panels but for a catalog with sample 
size comparable to the one used by Paraniape & Sheth 12012). The differences in the distribution from the two catalogs, as quantified by the p-values from the 
KS-test are indicated in each panel. 
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Fig. 5. — Distribution of p-values from KS-test carried out on the luminos- 
ity distribution of centrals (solid histogram) and satellites (dotted histogram) 
from the two catalogs car ried out on 1000 samples with a sample size com- 
parable to the one used by Paraniape & Sheth (2012). 

catalogs. 

It is well known that the luminosity of central galaxies 
depends upon the halo mass in which they reside. How- 
ever, it is also known that at the massive end, the lumi- 



nosity of central galaxies is a weak function of halo mass, 
e.g., based on two point statistics such as the projected 
galaxy-galaxy correla tion function, its dependence upon lu- 
minosity of galaxi es dZehavi et al.l 120051: IZheng et al.l 120071 : 



iZehavi et al.H201 ll) . and the projected galaxy-matter correla 
tion fu nction probed by the galaxy-galaxy lensing measure - 
ments (iMandelbaum et al.ll2006t ICacciato et alll2009l 120121) 
or other pro bes such as satellite kinematics ( More et al 
2009bLl2bl lb and subhalo abundance matchin g dMoster et al 



2010t iBehroozi et al.ll2010t lYang et al.ll2012h . This coupled 



with the fact that the satellite fraction is very low at the bright 
end, could be a reason why the differences in the magnitude 
distribution of the brightest galaxy between the two different 
catalogs are not that large. 

We note that this insensitivity of the magnitude distribu- 
tions to the un derlying halo occupation d istribution was also 
pointed out by lParanjape & Shethl (120121) . who suggested the 
use of two point statistics such as the luminosity-marked 
correlation function, in order to distinguish the two scenar- 
ios. Based on the radial dependence of the marked correla- 
tion function, they concluded that the galaxy luminosities in 
groups cannot be drawn from a global luminosity function. 
However, their analysis does not directly address whether this 
is due to a conditional luminosity function which varies with 
mass, or a result of environmental dependences of the lumi- 
nosity function. 



Magnitude Gap Statistics and the Conditional luminosity function 
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5. SUMMARY 

Recently, iHearin et al.l ( 120121) . suggested that the magni- 
tude gap between the two brightest galaxies in a given halo at 
fixed richness contains additional information about the halo 
mass. Their claim was based upon an analysis of the galaxy 
group catalog constructed from the SDSS by iBerlind et alj 
(2006). If correct, the magnitude gap information can be used 
to reduce the scatter in the mass-richness relation in galaxy 
clusters, which is important for the use of optically identi- 
fied galaxy clusters as cosmological probes. However, they 
cl aimed that their result is a t odds with the results presented 
in lParanjape & Shethj (120121) who investigated the distribution 
of magnitudes of the brightest a nd second brightest galax- 
ies, from the same group catalog. Paranjap e & Shethl (120121) 
showed that these magnitude distributions are consistent with 
the order statistics of the luminosities sampled from the over- 
all galaxy luminosity function independent of halo mass. This 
would imply that the magnitude gap just depends upon rich- 
ness and does not contain extra information about the halo 
mass. 

We have investigated both these studies within the frame- 
work of the conditional luminosity function (CLF), which de- 
scribes the halo occupation statistics of galaxies as a function 
of halo mass. The CLF and its variation with halo mass has 
been calibrated using observations of the abundance and clus- 
tering of galaxies, and the galaxy-galaxy lensing signal in the 
SDSS, and is consistent with results based upon the kinemat- 
ics of satellite galaxies and abundance matching. We have 
shown that if galaxies populate halos according to the CLF 
and if the luminosities of central and satellite galaxies are 
drawn from their corresponding CLF in an uncorrected man- 
ner, then the magnitude gap is expected to contain information 
about halo mass at fixed richness. We have presented analyt- 
ical expressions for predicting the magnitude gap distribution 
at fixed richness as a function of halo mass and verified these 
expressions using Monte-Carlo simulation of galaxy catalogs 
populated according to the CLF. 

We have shown that the magnitude distribution of the 
brightest and the second brightest galaxies show significant 
differences, between mock galaxy catalogs constructed by 
drawing galaxy luminosities according to the CLF and those 
constructed according to the luminosity function of galaxies. 
However, we have also shown that these differences cannot 
be meaningfully detected given the small sample size that 
iParanjape & Shethl (120121) use in their study. This shows that 



the magnitude distribution of the brightest and the second 
brightest galaxies is not the appropriate statistic to address 
the issue of how galaxies populate dark matter halos, at least 
given the current sample sizes. 

These results suggest that the apparent tension between 
the two studies is due to small sample size used by 
IParanjape & Shethl (120121) . The magnitude gap at fixed rich- 
ness can and does contain extra information about the mass 
of a halo. As the sample size of groups grows, even the lu- 
minosity distribution of the brightest and the second brightest 
galaxies will also be able to distinguish between the two sce- 
narios. In this paper, we have provided an analytical model 
based on the CLF to predict the magnitude gap distribution. 
We have also presented how the magnitude gap distribution 
can vary as some of our fiducial assumptions are changed. 

It is also important to note that the CLF of galaxies in clus- 
ters can also be directly observed, albeit as a function of opti - 
cal properties such as richness (see e.g.. lHansen et al.ll2009f) . 
Such observations, when combined with halo mass indicators 
such as weak lensing, can in turn be used to better constrain 
conditional luminosity function at the high mass end, which 
will help to constrain our model for the magnitude gap, at 
fixed richness. 

Finally, we would like to remark that we have assumed that 
the luminosities of the galaxies in every group are drawn from 
the conditional luminosity function in an uncorrected fash- 
ion. This assumption, however, needs to be thoroughly tested. 
For example, if bright satellite galaxies merge with the cen- 
tral galaxy due to dynamical friction, the central galaxy will 
become brighter at fixed halo mass, and the magnitude gap 
will correspondingly be larger. However, simultaneously the 
richness of the group will decrease (which will also cause the 
magnitude gap to be larger just due to the statistics of random 
draws from the CLF), thus making it difficult to disentangle 
the physical correlation from the effect due to changing rich- 
ness. 
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APPENDIX 

MAGNITUDE GAP WHEN SATELLITES ARE ALLOWED TO BE BRIGHTER THAN THEIR CENTRAL GALAXIES 

Based upon the analysis of offse ts of the line - of-sig ht velocities and projected position of the brightest galaxy relative to the 
mean of the other group members, Ski bba et al.1 (1201 ll) concluded that there is a significant chance that in a halo of given mass, 
a satellite galaxy is brighter than the central galaxy. In Fig. Q] presented in Section [3] we showed the magnitude gap distribution 
for the case when we allow satellite galaxies to be brighter than centrals. In this appendix, we provide analytical expressions for 
the magnitude gap (between the two brightest galaxies) distribution in this case. 

Note that since we allow for the possibility that the satellite galaxies can be brighter than the central galaxy, in a given halo, 
one of the following three mutually exclusive and collectively exhaustive cases may occur: (i) central galaxy still turns out to be 
the brightest, (ii) central galaxy turns out to be second brightest, (iii) central galaxy is not one of the two brightest galaxies. The 
probabilities corresponding to these three cases are given by 



Pi(M,N)= P C (L\M) [P s (< L|M)] dL, 



P 2 (M,N) = (N- 1) j P c (L\M)[l - P s (< L\M)] [P s (< L\M)] r 

P 3 {M, AO = 1 - P, (M, N) - P 2 (M, N) , 
respectively. The magnitude gap distribution can then be expressed as 

dL' P s (L'\M)P c (L'/f L \M) [P s (< L'\M)] (N - 2) 
+P 2 (M, N) (N - 1) dL' P S (L'\M) P c (L'f L \M) [P s (< L'f L \M)] 

dL' P S (L'\M) P s (L'/f L \M)[P s (< L'\M)] 

Anin 



dL, 
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(A3) 
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